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Abstract 

We improve tiie existing calculations of deep inelastic scattering to next 
to next to leading order in the following manner. First, we use the recently cal- 
culated values of the anomalous dimensions for moments with index n = 10, 12 
in ep scattering. Second, we use also recently calculated anomalous diemensions 
for n = 1,3, . . . , 13 for xF^ in i/N scattering, to extend the calculation to this 
process. Lastly, we use the determinantal constraints on the moments that follow 
from positivity of the structure functios to stabilize the fits. We find substantial 
improvement over previous results, getting the following numbers: 

A{nf = 4, 3 loop) =274 ± 20 MeV; '°°p(M|) = 0.1166 ± 0.0013; 
xVd.o.f. =92/(153 - 18) 
from ep and, from fN, 

Aijif = 4, 3 loop) =255 ± 72 MeV; al '°°p(M|) = 0.1153 ± 0.0063; 
xVd.o.f. =0.335/(64 - 7). 

The errors take into account experimental errors and higher twist effects among 
other estimated theoretical errors. 

We can also extrapolate the fit, made originally for the range 3.5 < < 
230 GeV^ to higher momenta, obtaining very good reproduction of experimental 
data up to ~ 5000 GeV^. ^From this, in particular, we can exclude super- 
symmetric partners (squarks or gluinos) for masses up to some 60 GeV. The fitted 
value of the coupling is now, including the data at high , 

'°°p(M|) = 0.1163 ± 0.0014 



Typeset with PnysMai^ 
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1 Introduction 

In two recent papers'^'^l (where we send for more details and and more extensive introduction and references) 

we have evaluated, and compared to experiment, the evolution of the structure function ^2(2;, Q^) in ep deep 
inelastic scattering (DIS) to next to next to leading order in QCD perturbation theory (NNLO): two loops in 
the Wilson coefficients, and three loops in the anomalous dimensions and strong coupling, as- We used the 
known values of the two loop Wilson coefhcients^^l and three loop anomalous dimensions'^^ for n = 2, 4, 6, 8 
for the moments for quark nonsinglet (NS), singlet (S) and gluon (G) functions, 

^0 



Mi(n;Q^)= / dxx''-^Fi{x,Q'); i = NS,S,G. (1.1) 
with 

F2{x, Q2) = Fs{x, Q^) + FNsix, Q'). (1.2) 

The evaluation used the method of Bernstein moments'^'^'. In this method one docs not compare 
with experiment the values of the ^{n, Q^) themselves; these cannot be evaluated reliably due to existence, 
for each value of Q"^, of only measurements of F2{x,Q'^) for a limited range of x. Instead, we compare 
theoretical predictions and experimental results for the Bernstein averages, 

Fnk{Q^)= f dxPnk{x)F2{x,Q^). (1.3) 

Jo 



Here the Pnk are the (modified) Bernstein polynomials, 

2r(n+|) 

"r(fc + i)r(n- fc + i) 



2{n-k)\r{n+l) (-1)' ^'-^^ 

r(k + \)r{n-k + i){^^i\{n-k-if ' 



These polynomials are positive and have a single maximum located at 

_ r(fc + i)r(n + |) 
r{k + \)r{n + 2y 

they are concentrated around this point, with a spread of 



\ n + 



r{k + i)r{n + i) 



r(fc + i)r(n + 2). 



2 



and they are normalized to unity, dxpnk{x) = 1. Therefore, the integral 



Jo 



nl 

dxpnk{x)F2{x,Q'^) 

Jo 

represents an average of the function F2{x) in the region Xnk — \^Xnk ~ a; ^ Xnk + \^Xnk'-, the values of 
the function outside this interval contribute little to the integral, as Pnk{x) decreases to zero very quickly 
there. So, by choosing suitably n, k we manage to adjust the region where the average is peaked to that in 
which we have experimental data. Finally, and using the binomial expansion in Eq. (1.3), it follows that the 
averages with the Pnk of F2 can be obtained in terms of its even moments: 

dxpr.kix)F2{x) + ^Z!(n-fc-0!^'™- 

This is then evaluated theoretically and the result compared to the experimentally evaluated averages. 

In the present note we improve on this as follows. First of all, for ep DIS, we incorporate the 
values of the three loop anomalous dimensions recently available'^l for n = 10, 12. Secondly, we also 
evaluate the evolution of the moments (Bernstein averages, in fact) for the structure function xF^ in vN DIS 
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(neutrino scattering on an isoscalar target). This we can do because the three loop anomalous dimensions 
for n = 1, . . . , 13 have also been evaluated rece ntlyl^l. So we have here a true NNLO calculation, as opposed 
to existing estimates where the value of the anomalous dimensions was extrapolated. 

The use of more moments means that we can also evaluate more Bernstein averages, and these more 
precisely, which in turn permits an important decrease of the statistical errors of the calculation. Because 
we have more moments, we can avoid the more suspect, lower data points, where uncontrollable higher 
twist terms or corrections due to the c quark mass (for example) are largest. So, for ep scattering, the energy 
range we choose is 3.5 GeV^ < < 230 GeV. 

We can also diminish the systematic errors by using the following improvement on the previous 
analysis. We profit from the fact that the structure functions are (proportional to) quark densities, hence 
positive-definite. This means that the moments must satisfy a set of very powerful determinantal inequalities, 
to be discussed in Sect. 2, which force a rigid behaviour of the moments and hence avoid spurious minima 
and spurious systematic deviations. This is particularly important for ep scattering, where we have eighteen 
input moments as free parameters. This is the second of the improvements that allow us a very precise 
determination of the strong coupling. We find. 



Q,(3 loop) (^2^ =0.1166 ± 0.0009 (statistical) ± 0.0010 (systematic) 
(for ep), 

and 

^(3 loop) (-^2^ =0.1153 ± 0.0041 (statistical only) 
(for uN). 

The results are compatible one with the other, and both with our "old" resultl^l for ep, 
Q,(3 loop) (^2-) ^0.1172 ± 0.0017 (statistical) ± 0.0017 (systematic) 
=0.1172 ±0.0024 {ep), 



(1.7) 

(1.8) 

(1.9) 



but the errors are now substantially smaller. 

The quality of the fit is so good that we can consider extending the analysis to higher values of Q^, 
where there are not sufficient experimental points to get Bernstein averages with only interpolation, so we 
have to resort to extrapolation to evaluate the integrals. This allows us to include data for ep scattering 
with momenta as high as ~ 5000 GeV^. We get two important results: i) If we include the new points 
in the fit, the value of obtained only changes to 

^(3 ioop)(^2 ^ ^ Q J163 ± 0.0010 (statistical error only); 



or we can simply extrapolate the theoretical fit obtained in the restricted range 3 - 230 GeV^. Then we 
find perfect agreement with the result obtained fitting the full range. These results show the stability and 
reliability of the calculations, as well as their suitability for high energy extrapolation. 

As a byproduct of this last analysis we can extend the results of ref. 2 to give bounds 



rriq, nig > 60 GeV 



for the masses of squarks or gluinos. 
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2 Bernstein Averages and Positivity Inequalities 



2.1 ep Scattering 

We consider as input parameters, to be fitted through the procedure to be described, the eighteen moments 
at a given reference momentum Qq, 

Hi{n, Ql), n = 2, 4, 6, 8, 10, 12; i = NS, S, G. (2.1) 

Actually, the momentum sum rule allows us to obtain iig{2,Q^) in terms of iis{2,Qf^), so we have only 
seventeen independent momenta; but, if we also fit the QCD coupling we have also A, so we still end up with 
18 parameters. We evolve the values of these moments with the QCD evolution equations (for the details, 
see ref. 1) so we get the moments to all values of Q'^: 

fj,i{n, Q^), n = 2, 4, 6, 8, 10, 12; i = NS, S, G. (2.2) 

With these momenta we form the theoretical Bernstein averages, which are then fitted to their experimental 
values. However, not any arbitrary input /Xi(n, Qq) are acceptable; they have to be such that, not only they 
fit the experimental data, but, at all Q^, they have to be the moments of positive-definite functions. 

It is possible to implement this restriction by imposing a set of determinantal inequalities, which 
are the necessary and sufficient conditions for the m{n, Q^) to be the moments of positive quark and gluon 
distribution functions. The details may be found in the Appendix, or in refs. 7,8. The inequalities may be 
written as 

A,[^i{n,Q^)]>0 1^ = 1,2,.... (2.2) 

where the are determinants formed with the moments. In our case, if we define Mk{Q^) = /x(2A; + 2, Q^), 
k ~ 0, . . . , 5, (for any of the NS, S, G cases) we can write (2.2) as the explicit inequalities that guarantee 
that the determinants 

/Ml M2 MgX 

det Ma M3 M4 , 
VM3 M4 Ms/ 

(Mo - Ml Ml - M2 M2 - M3 \ 
Ml - M2 M2 - Ms M3 - M4 I 
M2 - M3 M3 - M4 M4 - M5 / 

and all their principal minors are positive. Corresponding inequalities hold for the case when an odd number 
of moments is known. 

In principle, the inequalities (2.2) should be strictly satisfied, i.e., we should impose them on the fits. 

However, one should take into account that the calculation is not exact: higher twist, or JVNNLO effects are 
not included, and one should therefore allow a certain laxity. What we have done is to add to the chi-squared 
a penalty if some of the determinants are not positive. The penalty is obtained by multiplying the value of 
the negative determinants by 10^, and adding all these contributions to the chi-squared. This is done for all 
the values of Q^. In the actual fits, the sum total of the values of all the negative determinants is less than 
10~^, a quantity well within the tolerance due to the estimated errors stemming from higher order effects. 

2.2 vN Scattering 

The only difference in this case is that only odd moments are known now: 

fx{n,Q^)= [ dxx"-^xF3{x,Q^), n = 1, 3, 5, 7, 9, 11, 13. (2.3) 
Jo 

Therefore, to form the quantities Mi (see the Appendix) one has to consider F^ to be the structure function. 
The rest is as for the ep case, except that we have an odd number of moments now. 



- 3 - 



-J. SANTIAGO AND F. J. YNDURAIN- 



3 Computations and Results 



3.1 ep Scattering 

With the available values for the anomalous dimensions we can use the first six (even) moments. When 
evaluating the experimental input, wc have to calculate the average of the experimental F2 with Bernstein 
polynomials. In order not to rely on extrapolations, we follow the criterion stated in ref. 1 (where we send 
for details): we only include an average with pnk if we have experimental points covering the whole region, 
[xnk — ^Xnk,Xnk + ^Xnk]- For Bp Scattering this means that we can use the averages 

i^ir(O'), F^^^iQ'), F|o^p(Q2), 
F^^'{Q% F^^^{Q% F^^{Q% 

F^^^iQ'), F^^iQ^)- 

but of course not for each value of we have enough experimental information to evaluate all the moments. 
Altogether, we have 153 experimental averages.^ The values of the range from 3.5 GeV^ to 230 GeV^ 
The experimental input comes from ref. 9 for ep and ref. 10 for vN . The result of the fit may be seen in 
Fig. 1. 

Before presenting the numerical results, a few words have to be said both on the obtention of the 
experimental input and on the theoretical evaluation. As to the first, we remark that the calculation of an 
average such as (1.3), 

F„^^P(Q2)= f\xpnk{x)F;''''{x,Q% (3.2) 
Jo 

requires some interpolation and extrapolation of F^'^^'^^ {x, Q'^). To do so we have used two different methods. 
We separate F2 into a singlet and a nonsinglet part, writing 

F2 = Fs + Fns- 

In the first method we use, for each value of independently, a phenomenological expression for the F, 

Ff™■(:^)=(Ax-"■44 + C)(l-x)^ 
F^T'ix) =S.t"-5(1-x)^. 

We emphasize that, in the present paper, (3.3) is to be considered as only a convenient interpolation of the 
data to allow calculation of the integrals (3.2). In fact, the parameters A, B, C, v, ji are taken to be totally 
free, and assumed to be uncorrected for different values of Q^. Thus, no theoretical bias is induced in the 
dependence. 

The second method we consider is to use a polynomial parametrization of data. Then the integral 
involved in the Bernstein average is evaluated with the help of the parametric expression, for each value of 
Q^. A third method will also be used (for the nonsinglet piece); it will be discussed when we deal with vN 
scattering. As shown below, our results depend very little on the method of interpolation. 

We consider the first method to be the cleanest one; the results found using the polynomial 
parametrization are presented mostly to show the insensitivity of the evaluation to the method of obtaining 
the "experimental" averages f'^^'^\q'^). By varying F!^^"^ within its experimental error bars (including ex- 
perimental systematic errors) we obtain the averages, and their errors. The fit to these provides the values 
of our parameters and what we call our statistical error. 

The systematic error of our fit is really the resultant of the estimated theoretical errors, which we 
briefly discuss here: full details may be found in ref. 1. 



^ A point to take into account is that not all the averages are independent; because, for each we have only six 
moments, of tlie 12 possible averages only 6 arc independent. This can be taken into account by renormalizing the 
chi-squared scaling it to the number of independent experimental values; see again ref. 1 for details. This scaling 
has been taken into account, so the x^/d.o.f.presented here refers to the truly independent degrees of freedom. 
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Figure 1. NNLO fit to Bernstein averages of F2; TMC included. 



As stated before, we take as parameters of our fit the values of the moments at a reference momentum, 
/ii(n, Qq). For the fit we choose Qq = 8.75 GcV^. Wc may consider changing this to another value: the 
modification our results will suffer will be a measure of higher order effects. Other higher order effects are 
the JVNNLO corrections, that we estimate both as in ref. 1, and by using the four loop expression!^ ^1 for a^; 
higher twist corrections, that we incorporate phenomenologically by adding a contribution 

/.^^(n, Q2) = n{aA^/Q'')tXNs{n, Q^) (3.4) 

with a an unknown parameter, to be fitted, and expected to be of order unity. Next, we take into account 
target mass corrections (TMC), to order rn^/Q"^ by writing (sec again ref. 1 for details) 

(n, Q') = msin, Q') + Ifcil^/,;v5(n + 2, Q^). 

It should be noted that the TMCs for the singlet are negligible compared to those of the NS, which is why 

we only take into account the last. Wc remark that TMC arc incorporated into our best fit, but higher twist 
corrections or NNNLO ones are only used to estimate the theoretical uncertainty. 

Finally, we consider the matter of the number of quark flavours. The only quark threshold we cross 
is that of the h quark. Wc will be working with the moments ^i{n,Q'^), so we have only one momentum 
variable. We then split the range into the following two intervals: 

Q2 < ml (I); ml < (II). 
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Then, in region (I) we take n/ = 4 and in region (II), n/ = 5. The matching will be carried only for the 
coupling constant, following the standard prescription: to NNL0,[^^1 



/3o^+'log 



Po^)Lh + Snlo + Snnlo 



where 



= (&/ 

SnNLO = „n, j- 



an 

log 



.n/+l 



/?0 



,nf 



24 



Here, 



Lh = log [m\nf + 1)/A\nf)] , h = l3i/(3o 

and m(n/ + 1) is the pole mass of the (nj + l)th quark. For its vahic we take that of rcf. 13. 

It is important to realize that this method only provides an approximated evaluation of threshold 

To estimate the corresponding 



effects, which are somewhat distorted in the neighbourhood of — ml 

2 

mf 



6' 



i.e., to remove from the fit the 



error wc use an alternate method, which is to avoid the region 
corresponding experimental points. 

We present in Table 1 a compilation of the results obtained with our calculations at LO, NLO and 
NNLO, with TMCs taken into account; the fit to the data itself is shown, for the NNLO calculation (with 
TMC) in Fig. 1. Only statistical errors are shown in Table 1; systematic (theoretical) errors will be discussed 
below. We note that the errors given in Table 1 are "renormalized" to take into account the effective number 
of independent experimental points, as discussed before. 



Order 


A(nf = 4) 


«.(M|) 


xVd.o.f. 


LO 


226 ± 34 


0.131 ± 0.003 


130/(153 - 18) 


NLO 


279 ± 19 


0.1155 ±0.0014 


97/(153 - 18) 


NNLO 


274 ± 13 


0.1166 ±0.0009 


92/(153- 18) = 0.68 



Table 1 



The NLO corrections are very clearly seen in the fit: the x^/d-O.f. decreases substantially when 
including these. It also decreases when adding NNLO corrections, but the fit is so good already at LO that 
there is very little room for improvement. However the improvement is seen, particularly in that including 
NNLO corrections leads to a noticeable gain both in the quality of the determination of the coupling, 
and in the stability of the fits. Also the improvement in the quality of the fit with respect to that of 
refs. 1,2 is apparent; not only the errors are smaller, but the x^/d.o.f. has also substantially decreased, from 
79/(102 — 12) = 0.88 to the value 92/(153 — 18) = 0.68 we get now. In fact, the precision improves also to 
LO and NLO, thus emphasizing that the improvement, in particular the important gain in stability, is due, to 
a large extent, to the use of the positivity constraints that prevent the appearance of spurious minima. 

The estimated systematic errors, originating from various sources, are shown for the NNLO case in 
Table 2. It is interesting to remark that the use of the positivity constraints also improves the systematic, 
theoretical errors, which are now smaller than in the estimates of refs. 1,2, with the exception of the one due 
to estimated NNNLO effects. But this is because we also include here another source of error, namely, that 
due to cutting the expansion of as to three loops. We have. 
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Source of error 


Aim = 4; 3 loop) 


Ayl(n/ = 4;3 loop) 


AaJM^) 


No TMC 


279 


5 


0.0004 


Interpolation 


279 


5 


0.0004 


HT 


268 


6 


0.0004 


Quark mass effect 


269 


5 


0.0004 


Ql to 12 GeV^ 


279 


5 


0.0004 


AfNNLO 


265 


10 


0.0006 



Table 2 



Let us comment on the meaning of the different entries. No TMC means that we have not taken target mass 
corrections into account. The corresponding error is not included when evaluating the overall theoretical 
error because, since we take into account TMC in our central value, the error would be of order TMC^, 
or as X TMC, quite negligible. "Interpolation" means that we have used the polynomial interpolation to 
calculate the Bernstein averages. HT means that we have taken into account the presence of higher twist by 
adding a contribution like (3.4). The fitted value of the phenomenological parameter a is a = — 0.003±0.002. 
"Quark mass effect" mc^ans that we have cut off the b quark threshold, as discussed; the error in Table 2 
takes into account also the variations due to the error in the mi, mass, taken (with its value) from ref. 13. 
"Qq to 12 GeV^" means that we take the input moments to be defined at this value of the momentum, 
IJ,i{n, Qq = 12 GeV^), i = S, G, NS. Finally, JVNNLO means that we have fitted with the four loop formula 
for as- 

Composing quadratically systematic (theoretical) and statistical (experimental) errors we find the 
best result for the QCD coupling, 

A{nf = 4, 3 loop) =274 ± 13 (stat.) ± 15 (syst.) = 274 ± 20 MeV; 

a(3'°°P)(Mz) =0.1166 ±0.0009 (stat.) ± 0.0010 (syst.) = 0.1166 ± 0.0013. 

The corresponding central value for A(nf = 5, 3 loop) is of 193 MeV. 

It is to be noted that composing the "theoretical" errors as if they were independent leads to a certain 
amount of double-counting. Thus, the results should be independent of the value of Qq if the calculation was 
to all orders, so the two last errors in Table 2 are connected, and our estimate of the "systematic" (theoretic) 
error is really an overestimate. 

An important outcome of our calculations are the values of the moments of the quark and gluon 
densities, which should be useful as constraints for model-builders. These calculations we have performed in 
two ways. First, we can take directly the results of the fits. Secondly, we can impose with greater strength 
the positivity constraints at this initial value, = Qq by weighting the determinantal inequalities that 
ensure positivity more heavily there, so that the corresponding values of the moments arc more reliable. In 
both cases, the values of all the moments for the NS case, and of the lowest n (up to n = 8) for the quark 
singlet and gluon moments are essentially unchanged. Also, the values of A is left essentially unchanged, 
another proof of the stability of our results. The corresponding values of the moments are then, at our 
starting values of Q^ = Qo- 



Ql = 8.75 Ql = 12 

;xjvs(2) = 0.126 ±0.006 0.122 ±0.005 

;xjvs(4) = 0.0161 ± 0.0011 0.0150 ± 0.0010 

/"ivs(6) = 0.0036 ± 0.0005 0.0032 ± 0.0008 

UNsiS) = 0.0011 ± 0.0005 0.0009 ± 0.0004 

AtArs(lO) = 0.00041 ± 0.00017 0.00026 ± 0.00010 

Atjvs(12) = 0.00020 ± 0.00011 0.00011 ± 0.00006, 
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Ql = 8.75 








Ms(2) 


= 0.0521 ±0.0020 


0.0519 ±0.0015 






= (0.62 ±0.06) X 10-'' 


(0.77 ±0.06) X 


10^ 


-4 


MS (6) ^ 


= (0.24 ±0.08) X 10-'' 


(0.34 ±0.08) X 


10- 


-4 


Ms(8) : 


= (0.10 ±0.05) X 10-'' 


(0.23 ±0.07) X 


10- 


-4 




= (0.39 ±0.15) X 10-'"^ 


(0.10 ±0.05) X 


10- 


-4 


Ms(12) 


= (0.16 ±0.07) X 10-5 


(0.41 ±0.20) X 


10- 


-5 



and 

Ql = 8.75 Ql = 12 

Ijg{'2) = 0.208 ± 0.006 0.209 ± 0.004 

/xg(4) = (0.94 ± 0.008) x lO"-"^ (0.85 ± 0.08) x lO"-"^ 

^g(6) = (0.24 ± 0.06) X 10-3 (0.27 ± 0.05) x 10-^ 

yr^c(8) = (0.78 ± 0.28) X 10-^ (0.13 ± 0.03) x IQ-^ 

/xg(10) = (0.33 ± 0.12) X 10-" (0.50 ± 0.20) x 10"" 

AiG(12) = (0.10 ± 0.06) X 10-" (0.20 ± 0.10) x lO"" 

As discussed, the results are fairly stable; but we should emphasize that they do not include sys- 
tematic errors, which are expected to be largest for large n, especially for the quark singlet and gluon 
moments. 

If we do not overweight the reference moments we would have got the following values: 

Ql = 8.75 Ql = 12 

/XArs(2) = 0.126 ± 0.006 0.122 ±0.005 

UNsi'i) = 0.0160 ± 0.0010 0.0150 ± 0.0010 

I^Nsh) = 0.0036 ± 0.0005 0.0032 ± 0.0010 

Atjvs(8) = 0.0010 ± 0.0005 0.0009 ± 0.0004 

I^NsilO) = 0.00034 ± 0.00011 0.00026 ± 0.00009 

AtArs(12) = 0.00017 ± 0.00009 0.00012 ± 0.00007, 





Ql = 8.75 


02 = 12 






Ms(2) = 


0.0522 ±0.0015 


0.0519 ± 0.0008 




/^s(4) = 


(0.58 ±0.05) X 10-" 


(0.93 ±0.03) X 


10- 


-4 


Ms(6) = 


(0.27 ±0.08) X 10-" 


(0.54 ±0.02) X 


10- 


-4 


/^s(8) = 


(0.14 ±0.07) X 10-" 


(0.41 ±0.03) X 


10- 


-4 


Ms(lO) = 


: (0.41 ±0.12) X lO--"^ 


(0.94 ±0.10) X 


10- 


-5 


Ms(12) = 


: (0.22 ±0.10) X 10-5 


(0.55 ±0.09) X 


10- 


-5 



and 

Ql = 8.75 Ql = 12 

Mg(2) = 0.208 ±0.004 0.210 ±0.003 

/xg(4) = (0.86 ± 0.004) x IQ-^ (0.87 ± 0.04) x IQ-^ 

Mg(6) = (0.29 ± 0.06) X 10-^ (0.269 ± 0.013) x 10-^ 

^e(8) = (0.13 ± 0.08) x 10-^ (0.95 ± 0.11) x 10-" 

Mg(10) = (0.74 ± 0.08) x 10-" (0.53 ± 0.07) x IQ-" 

/Ug(12) = (0.24 ± 0.09) x 10-" (0.10 ± 0.05) x 10-" 

The difference with those reported above may be taken as a measure of the systematic errors involved in the 
calculations of the reference moments. 
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FlGURE 2. Fit to xFs with Bernstein averages. 



3.2 vN Scattering 

The treatment of xF^ in vN scattering is essentially similar to that of the nonsinglet part of F2 in ep 

scattering, except for the modifications indicated, which arc due to our having now only odd moments. We 
have now 64 Bernstein averages that can be calculated experimentally, and the range is from 7.9 to 125.9 
GeV^ 

We do not have now the momentum sum rule, but the Gross-Llewellyn Smith one, which we imple- 
ment. The results of the fit are shown graphically in Fig. 2. 

As for the numerical results, the estimated systematic errors are as for ep scattering; but, to test 

the dependence of the results on the paramctrizations used to get the experimental Bernstein averages, we 
have also calculated replacing the interpolating formula, that we now take as 

xFf"'-{x)=BxP{l-xY 

{p, jjL B freee parameters) by a cubic interpolation, requiring also the vanishing of the structure function at 
both a; = 0, 1. The results are given in the following Tables 3,4. 
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Order 


Aiuf = 4) 




xVd.o.f. 


LO 


217 ±78 


0.130 ± 0.006 


0.011 


NLO 


281 ± 57 


0.116 ± 0.004 


0.008 


NNLO 


255 ± 55 


0.1153 ±0.004 


0.38/(64 - 7) = 0.007 



Table 3 



The x^/d.o.f. is minute, the reason for which is obvious by looking at Fig. 2: clearly, the experimental 
errors are much overestimated. Nevertheless, it is satisfactory that the x^/d.o.f. of the fits decreases from 
LO to NLO, and even from this to NNLO. The systematic errors are shown for the NNLO case in Table 4: 



Source of error 


A{nf = 4;3 loop) 


AA{nf = 4;3 loop) 


Aa.(M|) 


No TMC 


298 


43 


0.0030 


Cubic Interpol. 


269 


14 


0.0010 


HT 


270 


15 


0.0010 


Quark mass effect 


240 


15 


0.0010 


to 12 GeV^ 


263 


8 


0.0005 


iVNNLO 


209 


46 


0.0037 



Table 4 



For the calculation with higher twist contributions, the value of the phenomenological parameter a 
is o = —0.14 ± 1.6, compatible (within the very large errors) with that found for ep. 
Composing errors as for ep we find, 

A{nf = 4, 3 loop) =255 ± 55 (stat.) ± 46 (syst.) = 255 ± 72 MeV; 

^(3 loop) (-^^^ =0.1153 ±0.0040 (stat.) ±0.0061 (syst.) = 0.1153 ± 0.0063. 

The values of the input moments are 

Ql = 8.75 Ql = 12 

)u^(l) = 2.69 ± 0.05 2.71 ±0.04 

/i,/(3) = 0.099 ± 0.018 0.095 ± 0.017 

(5) =0.016 ± 0.003 0.015 ±0.003 

/i^(7) = 0.0042 ± 0.00015 0.0039 ± 0.0013 

11^(9) = 0.0012 ± 0.0008 0.0011 ± 0.0005 

li^{n) = 0.00036 ± 0.00015 0.00032 ± 0.00011. 

M^(13) = 0.00015 ± 0.00011 0.00012 ± 0.00009. 

An important consistency check is the comparison of these with the nonsinglet ep moments (after interpo- 
lating to get from odd to even moments, and multiplying by the average charge factor, 5/18). We find, for 
Ql = 8.75, 



n 


{5/18)ii{,^N) 


ms{ep) 


2 


0.110 


0.126 ±0.006 


4 


0.010 


0.0161 ±0.0011 


6 


0.0020 


0.0036 ± 0.0005 


8 


0.0006 


0.0011 ±0.0005 


10 


0.00021 


0.00041 ± 0.00017 


12 


0.00009 


0.00020 ±0.00011 



and, for Qg = 12, 
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(GeV^) 
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0.1 








10 



1000 



100 
(GeV^) 

Figure 3. Extended fit to F2, to Q^=5000 MeV. 



n 






2 


0.109 


0.122 ±0.005 


4 


0.0097 


0.0150 ±0.0010 


6 


0.0019 


0.0032 ± 0.0008 


8 


0.00054 


0.00090 ± 0.0004 


10 


0.00020 


0.00026 ±0.00010 


12 


0.00005 


0.00009 ± 0.00004. 



The discrepancies are what one would expect from corrections due to perturbative, experimental and inter- 
polation effects, as well as the different isospin structure. 

3.3 Extrapolation to higher for ep Scattering. Bound on Squark and Gluino Masses 

As mentioned at the beginning of Sect. 3.1 we have, in the calculations carried till now, only considered fitting 
averages Fnk such that we have experimental points covering the whole interval [xnk — ^Xnk^Xnk + ^Xnk] 
where the corresponding Bernstein polynomial is concentrated; and this, to rely only on interpolations when 
evaluating the experimental averages. This is what limits the upper range of values of we can use 
to < 230 GeV^. Now, the results of our fits are so good that we may consider extending them to 
higher values of Q^. This will imply extrapolation to calculate the integrals, and hence larger errors for the 
experimental input, as well as a certain dependence on the parametrization we use to evaluate the Bernstein 
integrals. On the other hand, this will permit an important check of the possibility to extrapolate our results 
to higher Q^, an important issue in view of eventual applications to LHC studies. 
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To be precise, what we do is the following. We fit the recent HERA data on large Q"^ ep scattering 
for F2 with an expression like the one used before, 

FPi^«"-(a;) ={Ax-^-^^ + C){1 - x)\ 

We now admit a number of new averages, where the corresponding Bernstein polynomials arc peaked around 
X ~ 0. Specifically, we now admit averages provided only that at least six experimental points lie inside 
the interval [xnk — ^Xnk^Xnk + Ax„fe], (but we do not require that they extend to the edges) so that the 
fit be meaningful. In this way, we arc able to include experimental averages up to = 650 GcV^ for 
F^^^' , FgQ^' and up to = 5 000 GeV'^ for F^^' , F^^^'; other averages would get so large errors to make 
them meaningless. We then follow two strategies: i) We include the new points in the fit, which tests the 
stability of our results. If we do so we obtain the new results, 

A{nf = 4, 3 loop) =269 ± 14 MeV (stat.); 

Q,(3 loop) =0.1163 ± 0.0010 (stat.) ^^"'^^ 

with a x^/d.o.f. of /d.o.i. = 93.5/(174 - 18) = 0.60 slightly better than with the restricted fit.^ The 
agreement of (3.7) with (3.5) proves the noteworthy stability of our results. 

ii) Alternatively, we can extrapolate our results, with the central values as given in Eq. (3.5), plot 
the corresponding values for (the averages of) F2, and compare them both with data and with the curves 
obtained from the overall fit obtained following the procedure outlined in (i) above. 

The theoretical curves obtained with both methods, together with the corresponding experimental 
values, are shown in Fig. 3. The agreement of the curves found fitting also the new points, or simply 
extrapolating our previous results is so good you can hardly tell them apart in the figure. This shows very 
clearly the predictive power of our evaluation: from the input region, between 3.5 and 230 GeV^, to the 
5 000 GeV^ of Fig. 3. 

As an application of this, we improve on the bounds on squark and gluino masses obtained in ref. 2. 
Repeating the analysis carried there, we may conclude that the fact that the x^/d.o.f. actually decrease 
(slightly) when including high data implies that no new particles are excited in this range. This provides 
the bound 

rriq, rrig > 60 GeV 

for the masses of squarks or gluinos. 



4. Comparison with other calculations 

We next show a table comparing our results to previous determinations for as(M|). These are taken from 
the review by Bethke'^^l, except the Z — > jets one which is from ref. 15 and the one referred to as "our old 
result" which is of course that of ref. 1. We include only the processes where the theoretical calculation has 
been pushed to the NNLO level. 

^ We can even interpret the slightly smaller value of A obtained here as due to the greater weight that have now large 
data, so diminishing the contamination with higher twist effects. 
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Process 


Average 
or range [GeV]^ 


a.(M|) 


T^TQ- J/ Vt.'i 


1.58 


191+0.005 


"PlTC. ,, piT Q 
JJlO, i^, *jIjO 


3 


niio -U nmn 

[J.llZ zt U.UIU 


T decays 


(1.777)2 


n 1 1 Qi -U n nriQi 

U.iioi =t U.UUol 


e^e — > hadrons 


100 


1600 


U.lZo =t U.Uzo 


Z — » hadrons 


(91.2)2 


0.121+nni 

— U.UU4 


Z jets 


(91.2)2 


0.117 


Our result(:^Af; xF^) 


8 


120 


0.1153 ±0.0041 


Our oid result (ep) 


2.5 


230 


0.1172 ±0.0024 


Our result now (ep, I) 


3.5 


230 


0.1166 ±0.0013 


Our result now {ep, II) 


3.5 


5000 


0.1163 ±0.0014 



Table 5 



Here DIS means deep inelastic scattering, Bj stands for the Bjorken, and GLS for the Gross Llewellyn Smith 
sum rules. We have presented our two results here, the one for the smaller range (I, that we consider the 
more reliable) and the extended range one (II). The value given for Z hadrons assumes the Higgs mass 
constrained by 115 < Mh < 300 GeV, with the central value for Mh = 115 GeV. All determination are 
compatible with one another, within errors. 

We finish this paper with brief comments on other recent work on our subject. First, we have 
the evaluations of Katacv et al.'-'^^l who consider i/N scattering to NNLO. In their work, it is necessary to 
interpolate to obtain even moments and also to assume a Jacobi polynomial expansion for the structure 
function. They find the value A{nf = 4, 3 loop) = 332 ± 36, where some systematic errors are not taken 
into account. The result is compatible with ours, within errors. A second paper is that of van Necrvcn and 
Vogtl^^] who, however, are only interested in the consistency of the data with NNLO perturbation theory 
but, as far as we can see, no effort is made to extract the value of as- 

Appendix 

The constraints imposed on the moments of a positive function are based on the following result: 
Theorem (Lukacs). A polynomial positive in the interval [0, +1] may be written as 

P2n{x) = [ ^ ttkx'^j + X{1 - X) '^I'^j 
\fe=0 / \fe=0 / 

(even degree) 

P2n-i{x) = X l^akx''] +{l-x) ij2bkx''\ 

\fe=0 / \fe=0 / 

(odd degree). Sec e.g. Ahiezer and Krein, ref. 8, p. 35. 

For polynomials in x"^, we write P{x'^), P as before; likewise, xP{x'^) gives polynomials in odd powers 

of X. 

Positivity 
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Consider the moments ^ 

Mi= [ dif*(t), i = 0,l,...,N, (Al) 
Jo 

where the function ^ is positive {Hausdorff moment problem). We are going to describe necessary and 
sufficient conditions for the Mi to verify (Al). The conditions depend on whether N is even or odd, and are 
derived from the representation of polynomials positive in the interval [0, 1]. 

Theorem, i) Let N = 2m be even. Then, if Pjv is positive we write it as 

/ m \ / m— 1 \ 



vi=0 



^m— 1 \ /m—1 



ii) If = 2m - 1 is odd, 

PN{t)=t\^J2aifj +{l-t)\^J2bifj . (A2b) 

On integrating the P/v with \1/ we find the following conditions: 
i) For N even: 

m m—1 

jj^i — M^j^j^2] > 0; (A3a) 

i,j=0 i,j=Q 

for N =odd, 

m—1 m—1 

aittjMi+j+i + f}ibj[Mi+j - Mi+j+i] > 0. (A3b) 

i,j=0 i,j=0 

The problem is solved completely by using the following result: 
Theorem. The quadratic form 

n 
i,j=0 

is positive if, and only if, all determinants 



Al . 


• Afc \ 


A2 . 


• Afc+i 


fc+i 


A2fe / 



Afe[A] = det 



are positive for 2k < n. 
Connection with DIS 

For electroproduction, and F = Fns, Fs, Fg you know 

IJ.n= [ da;x"-^F(a;), n = 2,4,6,.... 
Jo 

Set then 

^ = a;^ n-2 = 2i, ^{t) = ^t-^^^F{t^/^)- 

then you have the Mi given by 

Mi= I dtt'^{t), i = 0,l,2... 
Jo 

For neutrino scattering, 

fin = dxx'^~^F3{x), n=l,3, 5, 

With n — 1 = 2i now, the rest as before, we get moments 



Jo 



(A4) 



Mi = / dxf^{t), i = 0,1,2... 
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